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We compute the spectral distribution of the quantum fluctuations of the vac- 
uum, amplified by inflation, after an arbitrary number of background tran- 
sitions. Using a graphic representation of the process we find that the final 
spectrum can be completely determined trough a synthetic set of working 
rules, and a list of simple algebraic computations. 
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The standard cosmological model U, in spite of the large number of successful predic- 
tions, cannot be extrapolated too far back in time. It is by now a consolitated opinion that 
our Universe, in its earlier epochs, should include a phase of inflationary (i.e. accelerated) 
evolution [fjj. Many different realizations of inflation are possible, however, and we have 
to face the problem of how to distinguish them through their possible phenomenological 
consequences. 

To this aim, we may note that the instability of the small fluctuations (of the metric 
and of the matter fields) is one of the main physical properties of a phase of accelerated 
evolution. As a consequence of this instability, small perturbations are amplified in a way 
which depends on the intensity and on the duration of inflation ||. Different inflationary 
models amplify perturbations with different spectral distribution, and the computation of the 
spectrum thus becomes an important tool to predict observable effects, and to discriminate 
eventually between various possible models of primordial evolution. 

The computation of the spectrum, performed accorting to the standard cosmological per- 
turbation theory 0, may become long and cumbersome even in the linear approximation, 
however, if the model of cosmological background is complicated. The aim of this paper is 
to present a set of working rules, based on a diagrammatic representation of the evolution 
of perturbations, allowing a quick estimate of the final amplitude and of the spectral distri- 
bution. The idea of a diagrammatic approach is not new |], but it was never implemented 
up to now in a complete and systematic way, so as to be generally valid for all classes of 
inflationary backgrounds. 

To be more precise, the diagrammatic method of computation that we shall propose here 
can be applied to all cosmological fluctuations ip described by the quadratic effective action 



where z is the external "pump field" , responsible for the amplification, and the prime denotes 
differentiation with respect to the conformal time 77. It is convenient, in this context, to 
introduce the canonical variable u = zip which diagonalizes the effective action, and which 



(for each mode k of the Fourier expansion, V 2 Wfc = —k 2 Uk)- We shall also assume that 
the background evolution can be separated into n + 1 different cosmological phases, with n 
transitions at 77 = r)i, i — 1, ...,n, and that in each phase the evolution of the pump fields 
(sufficiently far from the transition) can be parametrized in conformal time by an appropriate 
power ttj, namely Zj ~ \r]\ ai . It follows that, in each phase, the general solution of eq. (fj) 
can be written in terms of the first and second kind Hankel functions || as: 





satisfies the canonical evolution equation 
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4( V ) = \ V \ l > [A\(k)H^(x) + AL(k)Hff(x)] , m = \oi - 1/2|, x = k\ v \. (3) 

In order to fix our conventions, we shall label the n transitions in decreasing order for r\ 
ranging from — oo to +00, and we shall put the power on to the right (and Oj+i to the left) 
of rji, namely 

Zi+l ~ M ai+1 , V<Vi, Zi~\ri\ a \ 7]>r]i (4) 

It follows that the first (in order of time) cosmological phase is labelled by n + 1, the last 
one by 1. The normalization to an initial vacuum fluctuation spectrum |J thus imposes 
= 1, A™" 1-1 = 0, while the coefficients of the last phase A±(k), fixed by the continuity 
of ipk and ip' k at the transitions rji {not of Uk and u' k , see the discussion in ||), will determine 
the final spectral distribution of the amplified perturbations. The spectral energy density, 
in particular, is given by || dp/dlogk = (A;/ a) 4 (A;) | 2 /vr 2 or, in critical units and in terms 
of the proper frequency uj = k/a, 

The computation of the spectrum thus reduces, in general, to the problem of solving a 
linear non-omogeneous system of 2n equations, 

= ifitfa), = V£ 4 (*i), (6) 

for the 2n unknown quantities A±(k). The solution is straightforward, in principle; in 
practice, however, it is in general arduous to extract the relevant physical information from 
the exact solution, for any background with n > 2. In this paper we shall derive a set 
of prescriptions enabling an immediate estimate of the spectral coefficient A x _{k), through 
an approximate procedure which captures the essential features (amplitude and frequency 
dependence) of the spectrum, without solving explicitly the full system of equations. Such an 
estimate is based on the asymptotic expansion of the Hankel functions, and on the possible 
separation of the spectrum into n frequency bands, depending on the number of transitions 
which are truly effective for a given frequency mode. 

In order to introduce the basic ideas of our procedure, we shall start considering the 
simplest case of only one transition at 77 = 771. The normalized solution of eq. (0) is then 

ul = \ V \^H^(x), 



ul = \v\ 



A\HS(x) + A 1 _H^(x)} (7) 



(with [Mi 7^ /i 2 ), and the continuity of ipk at rji provides for the spectral coefficients the 
following exact solution H: 
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A - = f [" ( H '" H » - ) + {al - * 2) """" 

A \ = % h « >fl i?' - + («> - «o ^ < 

where the Hankel functions are evaluated ay xi = These coefficients satisfy the 

canonical normalization \A + \ 2 — \A_\ 2 = 1. 

The time scale rji defines the typical transition frequency, k\ = |?7i| _1 . For modes with 
k ^> k\ we can use the large argument limit of the Hankel functions, and we obtain A\_ ~ 0, 
A\ ~ 1: such modes are thus unaffected by the transition, modulo higher order corrections 
that are esponentially damped like e~ k l kx (see JF], for instance), and that we shall neglect 
in our approximation. For low frequency modes, k <C ki, we can use instead the small 
argument limit [[J, 

H^\x) ~ p^ + iq^x'^ - ir^x 2 '^ + s^x 2+ ^ + x -> 

Hf\x) ^p*^ -iq^x'^ + ir^x 2 -^ + s*^x 2+l1 + x-^0 (9) 

where the coefficients p,q, ... are complex numbers with modulo of order one (for later use, 
we have also included higher order corrections). The exact solution provides, in this 
limit, 

Al ~ C 1 1 xr" a " w + C\x^ +lxl + C\x^ - w + Clxl'^ 1 + ... (10) 

where 

= f (^2 - /il + «2 - «l) , 

Q2 = -\P»i%2 (-»2 -Hi-a 2 + ai) , 

C3 = ~f (^2 + /il - OS + «l) , 

c l = T [<?Mi r M2 ( 2 - y"2 + A*i - "2 + Qfi) + r w g w (-2 - p, 2 + /Ui - a 2 + (11) 

(and a similar expression for Ai). 

By recalling that /ij = |otj — 1/2 J , it follows that the first term of the expansion QTD|) is 
the leading one, for x\ — > 0. When = 0, however, we have to include the next-to-leading 
corrections. Taking into account all possible values of a%, a 2 , and truncating the expansion 
(PHI) to the lowest order term with non-vanishing coefficients, we find that for ati 7^ «2 there 
are four different possibilities, corresponding to four different spectral amplitudes: 



-Mi 

■ 1 xr Ml+1 , (12) 



ct 2 > 1/2 or «i > 1/2, 


A\ 


~ C\x^ 2 


«i > —1/2, «i > «2, 


A l _ 




a 2 > -1/2, a 2 > Qti, 


Ai 


~ Cjxf a: 


a 2 < -1/2, a x < -1/2, 


Ai 
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in agreement with the results first obtained in § . 

If one of the two powers «i, a 2 equals 1/2, the corresponding Bessel index is fi = 0, and 
the small argument expansion (|9]) is to be replaced by 



Hq 1 ^ (x) ~ po + iqo log x — ir x 2 log x + s^x 2 + . . . , 
Hq^ (x) ~ p* Q — iqo log x + ir x 2 log x + s* Q x 2 + . . . , 



x 
x 



0. 
0. 



Suppose, for instance, that a.\ = 1/2. The exact solution (H) is now approximated by 



A 



C\x^ 2 \ogx x + C^ 2 + 



where 



7T 
1 



Wo — M2 



1 

a 2 H — 
2 



(13) 



(14) 



(15) 



As C 1 is always nonzero, the leading term in eq. (|T4]) is the firs one for a>2 > 1/2, and the 
second one for ai < 1/2. We may thus include also the value a = 1/2 in the general rules 
(|j~2l), provided we take into account the prescription 



X \(J,=0 



logs, 



x \n=0 



1. 



(16) 



The above computation for a single background transition can be easily iterated for a 
cosmological model containing two or more transitions. We have solved the general case with 
n transitions, and we have found the remarkable result that the vanishing of the leading term 
of the asymptotic expansion - and then the particular spectral behaviour of the solution - 
depends only on the kinematic powers «j of two phases: the one preceding the first transition, 
and the one following the last transition. Such a result is in agreement with the well known 
phenomenon of "freezing" of perturbations 0, and with the general duality properties |§ of 
the action ([[]). 

Using the above result, it becomes possible to write down a recurrent expansion for the 
spectrum after n transition. To this aim, it is convenient to represent the whole amplifi- 
cation process with a simple diagram, in which we insert a vertical line in correspondence 
of each transition. The height of the i-th line is proportional to the associated transition 



frequency, ki 



\Vi 



It become possible, in this way, to identify at a glance the various 



frequency bands (k) of the spectrum, according to the number of transitions (with ki > k) 
from which a given band is significatively affected. Note that, for rj growing from minus in- 
finity, the height of the vertical lines may grow monothonically up to a maximum transition 
frequency (corresponding to a minimum time scale), and is then monothonically decreasing 
for n running towards plus infinity. 
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FIG. 1. The n frequency bands of the spectrum for a background with n transition frequencies, 
arranged in growing order. 



Let us start with the case in which the transition frequencies ki = are arranged in 

growing order from the left to the right (see Fig. 1). The frequency band with k > k\ is not 
amplified and we shall disregard it. The band number 1 of the diagram [k 2 < k < k±) will 
be affected by one transition, the band number 2 (k 3 < k < k 2 ) by two transitions, and so 
on. The iteration of the matching procedure used for a single transition leads to a recurrent 
expression for the coefficient A_, where the ^-dependence is fixed by the first and last phase, 
and the amplitude is fixed by the continuity at the transitions. We can write, in particular, 

\A l _(k)\~xl 21 xT\ k 2 <k<h 

„731 _713 



' r , \ 721 / r ,\ 712-713 
xTxri-) (-) , h<k<k 



x 2 J \x 2 J 



/ T 1 \ T™ 1 / Tl \ 7ln-7l,n+l 

7«+l,l-r7 1 >»+ 1 ( n ~ l ) ( Zl \ 

X n / \x n/ 



2 



— X n X l 



ftiyvfiy^way-fE.^. k<kn (17) 

Xi J \XiJ \X 2 J \X 2 J 

(modulo a numerical factor of order one, that can be computed from the exact solution, but 
that we shall neglect for our purpose of a quick approximate estimate). 

The powers 7^ of eq. Ql7|) depend (in ordered way) on «j, according to the same rules 
of eq. (0), with the only difference that when «j and are not contiguous (i.e., «j 7^ a^+i), 



they can also have the same value. In that case /ij = and the lowest order non- vanishing 
term of the expansion leads to x~^{Xf.* if < x^, or to x^x^ 1 if Xi-\ > Xk (see eq. (|TTD, 
where C\ = — C\ 7^ for fii = fi 2 ). By summarizing the working rules for computing the 
powers 7^ we must then distinguish two cases. In the case a« = we have: 
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FIG. 2. The n frequency bands of the spectrum for a background with n transitions frequencies, 
arranged in decreasing order. 



OLi = a k , 



e ik = sign log 



Ik 

Xj-l 



sign log 



Otherwise (a* 7^ a*.) we have, according to eq. (|T2|): 



«i > 1/2 or « fe > 1/2, 
«i > -1/2, a, > a k , 
a k > -1/2, a k > a h 
on < -1/2, a k < -1/2, 



7ifc l^i ; 

7ifc = — A*i, 
7ifc = — A** + 1, 



7fci A^fcj 

7fcj = -[i-k + 1, 



(19) 



Note that j ik 7^ 7^, but that 7^(0;,, a^) = 7ifc(afc, otj). Note also that, in the limiting case 
in which //£ = or ri k = 0, we can take into account the logarithmic corrections according 
to the prescription (IB). However, as they are usually negligible in realistic models, we shall 
neglect the log corrections in our first estimate of the spectrum, using the simple rule: 



lik \ ^,a k 



0. 







1 


Iki (2'°*) 


= -fik = 


a k -~ 



(20) 



for any k. 

At this point, one remark is in order. The spectral distribution determined through 
the above procedure applies - by construction - only sufficiently far from the transition 
frequencies, k -C k{. Near the transition the spectral slope may change with respect to 
the asymptotic regime, but always in such a way as to guarantee the continuity of the 
spectrum at k — k { . In eq. (|17D, however, the asymptotic expression has been extrapolated 
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up to hi and, as a consequence, its amplitude is normalized by continuity. This is certainly 
allowed - within our approximations - when the leading term of the expansion ( [TO] ) is 
nonzero. When a, and ol\ are both smaller than 1/2, and the leading term is vanishing, this 
extrapolation may introduce an error in the asymptotic amplitude of the spectrum, which 
is however of order one - and thus compatible with the degree of accuracy required for our 
present estimate - provided the amplification of the corresponding band is not too large, 
i.e. provided F{ki) = {k i /ki^- i ) lil {ki/ki) lu ^ 1. If F 1, on the contrary, the correct 
asymptotic amplitude of the branches with k < ki is obtained by renormalizing the results 
of eq. (0) by the factor F" 2 ^). 

A recurrent expression, similar to eq. (|17D , is also obtained for the complementary 
situation in which the n transition frequencies ki are arranged in decreasing order, for rj 
ranging from minus to plus infinity (see Fig. 2). There are still n frequency bands, defined 
by < k < ki, i = 1, ...,n, and the spectral coefficients can be approximated as follows: 

~ xl n+1 * n xl n - n+1 , k n _ x < k < k n 

k-n—2 ^ k < k n —i 





7n,n+l / \ 7n+l,n— 7n + l,n-l 

/ Xn \ 



\-En- 1 / 

7„+l,„-l-7„+l,n-2 /£ 2 \72,n+l / x \ 7n+l,2-7n+l,l 
X n - 2 J \X n -2j "xXtJ \XiJ 

k < ki (21) 

where 7^ are given again by eqs. (|T8|), (p!9|). Again, when the leading term is vanishing, 
the asymptotic amplitudes are correct provided F(ki) = (/i;i//i;i+i) 7i+1 '™ +1 (A;i/A; n ) 7n+1 '* +1 < 1, 
otherwise they are to be renormalized by the factor F 2 (ki). 

The above results can be summarized by a set of prescriptions, allowing an automatic 
computation of the spectrum, once the relevant diagram is plotted. For a synthetic formula- 
tion of such prescriptions it is convenient to define, for each "leg" of the diagram representing 
a transition, the so-called phase of exit and phase of re-enter. More precisely, for any leg rji 
placed to the left of the highest one, we shall define phase of re-enter the one directly to the 
right of the last leg crossed by the transition frequency ki = For any leg rjj placed to 

the right of the highest one, we shall define phase of exit the one directly to the left of the 
first leg crossed by the transition frequency kj = (see Fig. 3). 

With the above definitions, the diagrammatic rules for the computation of the spectrum 
can now be synthetized as follows. 

• We plot the diagram for the given model of background 
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FIG. 3. A graphic representation of the exit phase for the "leg" rjj, and of the re-enter phase 
for the "leg" rji. 



We choose the frequency band we want to compute, and we single out the relevant 
transitions. 

For any relevant leg rji placed to the left of the highest one we insert the amplitude 
factor: 

-?±) (^) . (22) 

For any relevant leg rji placed to the right of the highest one we insert the factor: 

ft ■ i i \ 7i+l,e f T 1 \7e,i+l— lei 

(t 1 ) (¥) ' (23) 

where the labels "e" and "r" denote, respectively, the exit and re-enter phase (note 
that, according to these rules, the highest leg does not contribute to the amplitude). 



We add the overal factor determing the frequency dependence of the spectrum 

j f x e i 



xj +1 ' e xj ej+1 , (24) 



where the labels "/" and denote, respectively, the first and last (from the left) 
relevant legs for the band we are considering. 

• We compute, finally, the powers 7^ according to eqs. (0), (|19D and, if needed, we 
renormalize the amplitudes through the factors F~ 2 or F 2 , as discussed before. 

As a simple application of such a method, we shall compute here the spectrum for the 
amplification process represented by the diagram of Fig. 4, by assuming for the various 
phases the following particular values: 



9 



a 4 


a 3 


a 2 


• 

a! 


113 112 


ill 



FIG. 4. An example of diagram representing the amplification of tensor metric perturbations, 
in a non-minimal model of string cosmology. 



a x = 1, a 2 = 1/2, a 3 = 2, a 4 = 1/2. (25) 

With such a choice of the kinematical powers, the diagram of Fig. 4 may represent (in 
the Einstein frame) the amplification of tensor metric fluctuations in a non-minimal model 
of string cosmology @||, which includes an initial dilaton-driven phase (rj < 773), a first 
intermediate, high-curvature string phase (773 < 77 < rj 2 ), a second intermediate dilaton 
phase with decreasing curvature (772 < 77 < rji), and a final radiation-dominated phase, with 
constant dilaton (77 > rji). In the limit 772 — > 77! the model reduces to the minimal one, in 
the limit 773 — > 772 one recovers instead the non-minimal model discussed in ||. 

Let us consider, for instance, the (lowest frequency) band number 3 of Fig. 4, corre- 
sponding to a "three-leg" transition (for this band, obviously, 773 is the first leg and 77! is 
the last one). For the leg 773 the re-enter phase is a 2 , for the leg 771 the exit phase is a 4 . 
Following the rules listed above we obtain the spectral coefficient 

The computation of the powers 7^, according to eqs. (pT8|), (|19D, (|20D, gives: 

732(2,!) = — 723(2,!) =0, 

724 (|, I) = 0, 742 Q,I) =0 , 

74i(il) = 0, 714 Q>l)=-^> (27) 

where the ordered labels of 7 refers to the various phases, and the numbers enclosed in round 
brackets refer to the particular numerical values of the corresponding kinematical powers. 
The spectral coefficient (|26|) thus becomes 
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) 



3/2 



k 



) 



1/2 



(28) 



i 



and the associated energy distribution (|5|), 




3 



) 



3 



fi(w,t) 



3tt ikf 2 // 2 



(29) 



reproduces (modulo logarithmic corrections) the well known cubic slope associated to the 
dilaton-radiation transiton, and is in agreement with the results of [(§] and [[J. With the 
same procedure we can easily estimate the spectrum for the other bands of Fig. 4. 

In summary, we have shown in this paper how to obtain a quick estimate of the cosmo- 
logical spectra using a simple method, based on a set of effective prescriptions, and on a 
diagrammatic representation of the amplification of perturbations. We believe that such a 
method does not represent only a mathematical "curiosity", but may have useful applications 
to the study of various inflationary scenarios, that will be discussed in future papers. 
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